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Abstract 

The gluon dominance model offers a description of multiparticle production in 

e + e~ - annihilation and proton-proton collisions. The multiplicity distributions in 

{Sj ■ e + e~ annihilation are well described. The energy dependence of model parameters 

gives the dynamic parton stage and hadronisation picture. It is shown that this 

^ ', model has confirmed oscillations in sign of the ratio of factorial cumulant moments 

over factorial moments of the increasing order. 

The collective behavior of secondary particles in pp-interactions at 70 GeV/c is 

studied in the project "Thermalization". An active role of gluons is shown in the 

multiparticle dynamics. This paper gives a simple thermodynamic interpretation of 

interactions mentioned above. 
(N ■ 
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§ ; 1 Introduction 

in ■ 

At present to investigate and construct a contemporary picture of nuclear matter structure 
requires to develop new methods and approaches. 

Since 80's the Quark-Gluon Plasma conception has undergone a lot of changes after 
the experiments carried out at CERN (SPS) and BNL (RHIC). Different approaches are 
used to explain extraordinary phenomena in behavior of the new matter produced at high 
>• . energy of nuclear collisions [I]. 

Still there is no single theory nor the model that could explain all the results obtained 
at RHIC and SPS. We would like to find a solution for this difficult problem by using 
multiparticle production (MP) in hadron and nucleus interactions. Up to now the nature 
of " soft" hadronic events has not been fully understood. 

A new way to investigate MP at high energy is offered in this work by means of 
construction a model based on the multiplicity distribution (MD) description using the 
QCD and phenomenological scheme of hadronisation. The model description of MD in 
e + e~ annihilation is given in section 2. The application of this model approach to pp- 
interaction can be found in section 3. 
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2 MD for e + e -annihilation at high energies 

The e + e~ annihilation is one of the most suitable to study MP. It can be realized through 
the formation of virtual 7 or Z°-boson which then decays into two quarks: 

e+e" - (Z°/ 7 ) - qq- (1) 

1 Talk given at the XVII International Baldin Seminar " Relativistic Nuclear Physics and Quantum 
Chromodynamics" . JINR. September 27 to October 2, 2004, Dubna, Russia. 
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The e + e~-reaction is simple for analysis, as the produced state is pure qq. It is usually 
difficult to determine the quark species on the event-by-event basis. The experimental 
results are averaged over the quark type. Perturbative QCD (pQCD) may be applied to 
describe the process of parton fission (quarks and gluons) at big virtuality, because strong 
coupling a s is small. This stage can be called as the stage of cascade. When partons get 
small virtuality, they change into hadrons, which we observe. At this stage we can not 
apply pQCD. Therefore phenomenological models are used to describe hadronisation in 
this case. 

Parton spectra the quark and gluon fission in QCD were studied by K. Konishi, 
A. Ukawa and G. Veneciano. The probabilistic nature of the problem has been estab- 
lished [2] while working at the leading logarithm approximation and avoiding IR diver- 
gences by considering finite x. Studying MP at high energy we used ideas of A. Giovan- 
nini [B] to describe quark-gluon jets as Markov branching processes. Giovannini proposed 
to interpret the natural QCD evolution parameter Y = ^ln[l + a s 61n(^-)] , where 
2nb = |(lliV<7 — 2Nf) for a theory with N c colours, Nf flavours and virtuality Q as the 
thickness of the jets and their development as the Markov process. 

Three elementary processes contribute into QCD jets: (1) gluon fission; (2) quark 
bremsstrahlung and (3) quark pair production. Let AAY be the probability that gluon 
will convert into two gluons in the infinitesimal interval AY, AAY be the probability that 
quark will radiate a gluon, and BAY be the probability that a quark- ant iquark pair will 
be produced from a gluon. A. Giovannini constructed a system of differential equations 
for generating functions (GF) of quark and gluon jets 

^ = AQ^(Q^ - 1), ^ = A(Q^ 2 - QW) + B{QM - Q {9) ) 

and obtained explicit solutions of MD in the case B = (process of quark pair production 
is absent) 

P P (Y) = ( ^ p \ P P P (Y) — kp(^P + !)••• {hp + m— 1) f m \ ( k p ^ 



k p + m J "' ml \m + kpj \k p + mj 

(2) 

where k p = A/A, m = k p (e AY — 1) is the mean gluon multiplicity. These MD are known 
as negative binomial distribution (NBD). The GF for them is 

oo 

Q®{z,Y) = £ z m P m (Y) = [l + m/k p (l - z)]- hp . (3) 

m=0 

Two Stage Model (TSM) |lj was taken (J2J) to describe the cascade stage and added 
with a sub narrow binomial distribution for the hadronisation stage. We have chosen 
it basing on the analysis of experimental data in e + e~- annihilation lower than 9 GeV. 
Second correlation moments were negative at this energy. The choice of such distributions 
was the only one that could describe the experiment. We suppose that the hypothesis of 
soft decoloration is right. Therefore we add the hadronisation stage to the parton stage 
for the sake of its factorization. MD in this process can be written as follows: 

Pn{s)=Y.Pn l Pn{m, S ), (4) 
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where P^ is MD for partons (|2*|). P^(m, s) - MD for hadrons produced from m partons at 
the stage of hadronisation. Further we substitute variable Y on a center of masses energy 
yfs. MD of hadrons P^ formed from one parton and their GF Q?(z) are [3] 
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where C^f p - binomial coefficient, and N p (p = q,g) have a sense of mean multiplicity 
and maximum of secondary hadrons are formed from parton at the stage of hadronisation. 

MD of hadrons in e + e~ annihilation are determined by convolution of two stages 
(cascade and hadronisation) 



Pn(s) = £ P , 
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where 2 + m is the total number of partons (two quarks and m gluons). 

We introduce parameter a = N g /N q to distinguish the hadrons, produced from quark 
or gluon. Also we have carried out simplification for designation iV = N q , n h = Ti h q . 
Introducing expressions (J2J), © in © and differentiating on z, we obtain MD of hadrons 
in the process of e + e~ annihilation in TSM 
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The results of comparison of expression (J7J) with experimental data [3] are shown in 
Figs. 1-2. We have obtained that MD in TSM (solid curve) describe well the experimental 
e + e~-data from 14 to 189 GeV . The mean gluon multiplicity m has a tendency to rise, 
but lower than the logarithmic curve. Values k p remain ~ 10 at almost all energies. One 
of the most interesting physical senses of this parameter is temperature T[7]: T ~ k~ x . 

The next picture of the hadronisation stage is discovered in conformity with parameters 
of the second stage: N, n h and a . The first parameter N determines the maximum 
number of hadrons, which can be formed from quark while its passing through this stage. 
We can imagine that fission is continuous but process (3) (formation qq pair) becomes 
comparable with the other ones (1),(2). We can not reveal a steady energy rise or fall for 
N. 

The second parameter n h has a sense of the mean hadron multiplicity from quark at 
the second stage. We have found out the tendency to a weak rise with big scattering. 
The value of n h is about 5 — 6 in the research region. A possible explanation of these 
rocks for iV and n h : only two initial quarks exist among a lot of newly born gluons at the 
beginning of hadronisation. 

The last parameter a was introduced to compare the quark and gluon hadronisation. 
It is equal to 0.2 with some deviations. If we know a, then we can determine N g = aN 
for gluon. It is surprising that gluon parameters remain constant without 
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considerable deviations in spite of the indirect finding: N g ~ 3 
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1 (Fig. 3-4). 



Therefore we can confirm the universality of gluon hadronisation. The fact that a < 1 
shows that hadronisation of gluons is softer than that of quarks. 
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It was shown jH] that the ratio of factorial cumulative moments K q over factorial 
moments F q changes the sign as a function of the order. We use MD formed in TSM to 
explain this phenomenon. F q and K q are obtained from the relations 

oo q—l 

F q = n(n - 1) . . . (n - q + l)P n , K q = F q ~Y / C^K^F,. (8) 

n=q i=\ 

The ratio of their quantities is H q = K q /F q . The generating function for MD of hadrons 
(J7J) in e + e~ annihilation Q(z) is the convolution 

QW = E P 9 M (zTQfc) = Q 9 (Qf(z))Q 2 q (z). (9) 

m=0 

We calculate F q and K q in TSM, by using and the sought-for expression for H q will 
be U 

£ k p am(am - i) ... (am - ^)(=^) m i - 2(-l)^ 

jt _ m=l (10) 

q E (2 + am) (2 + am - ±) . . . (2 + am - ^-)P m ' 

m=0 

The comparison with experimental data jH] has shown that ()10|) describes the ratio of 
factorial moments (Fig. 5). The minimum is seen at q = 5. We have obtained that in the 
region before Z° (91.4 GeV), H q oscillates in the sign only with the period equal to 2 and 
changes the sign with parity q. At higher energies the period is increased to 4. It can be 
explained by influence of a more developed cascade of partons with narrow hadronisation. 



3 MD in ^-interactions 

The study of MD in pp interactions is implemented in the framework of the project 
"Thermalization". This project is aimed at studying the collective behavior of secondary 
particles in proton-proton interactions at 70 GeV/c 9J. On the basis of the present 
understanding of hadron physics, protons consist of quarks and gluons. After the inelastic 
collision the part of the energy of the initial motive protons are transformed to the inside 
energy. Several quarks and gluons become free. Our model study has shown that quark 
branching of initial protons in pp interactions is almost absent from 70 to 800 GeV/c. MP 
are realized by active gluons. Domination of gluons was first proposed by S. Pokorski 
and L. Van Hove |10j . 

Our choice of the MP model is based on comparison with the experimental partial 
cross section u(n ch ) in pp interaction at 70 GeV/c on the U-70 accelerator [II] and the 
present picture of strong interactions. 

At the beginning of 90s a successful description of MD at this energy was realized 
by the quark model (Fig. 6) [HJ. This model suggests that one proton quark pair, two 
pairs or three can collide and fragment into hadron jets. MD in quark jets were described 
by Poisson. Second correlation moments of charged particles for MD in this model will 
be always negative. It is known they are become positive at higher energies. In this 
model gluons are absent. The calculation by the MC PHYTHIA code has shown that 
the standard generator predicts a value of the cross section which is in a reasonably good 
agreement with the experimental data at small multiplicity (n c h < 10) but it underesti- 
mates the value a{n ch ) by two orders of the magnitude at n ch = 18 (Fig. 6). 
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We have managed to build a scheme of hadron interactions to describe MD with the 
quark-gluon language as well as to investigate the high multiplicity region. The mentioned 
models are very much sensitive in this region. 

We consider that at the early stage of pp interactions the initial quarks and gluons 
take part in the formation of quark-gluon system (QGS). They can give branches. We 
offer two model schemes. In the first scheme we study hadroproduction with account of 
the parton fission inside the QGS and build the two stage model with branch (TSMB). 
If we are not interested in what is going inside QGS, we come to the thermodynamical 
model (TSTM). Onward we name models involving active gluons into hadroproduction 
as the gluon dominance models (GDM) [TB*] . 

We begin our MD analysis with the scheme of branch. MD for quark and gluon jets 
may be described NBD and Farry distributions jH], accordingly. On the hadronisation 
stage we have taken a binomial distribution As in TSM we have used a hypothesis 
of soft decoloration for quarks and gluons at their while passing of this stage and add the 
hadronisation stage to the branch one by means of factorization 



where P n (s) - resulting MD of hadrons, P p - MD of partons (quarks and gluons), P^{m) - 
MD of hadrons (second stage) from m partons. Generating function (GF) for MD in 
hadron interactions is determined by convolution of two stages: 



where Q H and Q p are GF for MD at hadronisation stage and in QGS. 

At the beginning of research we took model where some of quarks and gluons from 
protons participate in the production of hadrons. Parameters of that model had val- 
ues which differed very much from parameters obtained in e + e~- annihilation, especially 
hadronisation parameters. It was one of the main reasons to refuse the scheme with active 
quarks. After that we chose the model where quarks of protons did not take part in the 
production of hadrons, but remained inside of the leading particles. All of the newly born 
hadrons were formed by gluons. We name these gluons active. They could give a branch 
before hadronisation. 

It is important to know how much active gluons are into QGS at the first time after 
the impact of protons. We can assume that their number may change from zero and 
higher. It is analogous with the impact parameter for nucleus. Only in the case of elastic 
scattering the active gluons are absent. The simplest MD to describe the active gluons 
formed in the moment of impact is the Poisson distribution Pk = e~ k fk\, where k and 
k are the number and mean multiplicities of active gluons, correspondingly. 

To describe MD in the gluon cascade formed by the branch process of k active gluons, 
we have used the Farry distribution [3] 
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Q(s, z) = £ PZW {Q h (z)T = Q P (s, Q H { z ))i 
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at k > 1 - (|T3"j) and at k = 1 - ()14|) . m and m are the number of secondary gluons 
and mean multiplicites of them (averaged to all gluons). Expressions (I13)) - (jl4)) have been 
obtained from the assumption about the independent branch of gluons 
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In the case k = (the impact was elastic and active gluons are absent) MD of hadrons in 
pp-scattering is equal to -P 2 e '( s ) = e ~ k ■ 

On the second stage some of active gluons may leave QGS and transform to real 
hadrons. We can name that gluons "evaporated". Let us introduce parameter 5 as the 
ratio of evaporated gluons, leaving QGS, to all active gluons, which may transform to 
hadrons. Our binomial distributions for MD of hadrons from the evaporated gluons on 
the stage of hadronisation are 
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__ h \ SmN- (n-2) 



n \ 



(16) 



In this expression the gluon parameters are n h and iV (without index "g") which have 
the same meaning that of the quark parameters. An effect of two leading protons is also 
taken into account. GF for MD (fTp^l has the following form: 



Qm = (0? 



Sin 
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(17) 



MD of hadrons in the process of proton-proton scattering in two stage gluon model 
with branch (TSMB) is 
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In comparison with experimental data [TJ the numbers of gluons in sums on k and m 
were restricted by values MK and MG as the maximal possible number of gluons on the 
transition. For comparison we have taken the data at 69 GeV/c because they do not differ 
from data at 70 GeV/c [TTj j. x 2 /ndf in are equal to about ~ 1/3 at 70 GeV/c and ~ 5 
at 69 GeV/c and the parameters are similar. We obtained iV = 40(fix), m = 2.61 ± 0.08, 
5 = 0.47 ± 0.01, k = 2.53 ± 0.05, n h = 2.50 ± 0.29 from the comparison with [TT] • We 
can conclude that the branch processes are absent, since parameters m and k are equal 
to the errors. The fraction of the evaporated gluons is equal to 0.47. A maximal possible 
number of hadrons from the gluon looks very much like the number of partons in the glob 
of cold quark-gluon plasma of L.Van Hove ^3]. The gluon branch should be very active 
inside of QGS. At the fixed parameter of hadronisation n h equal to 1.63 (see below the 
thermodynamic model) the fraction of the evaporated gluons 5 is about 0.73. After the 
evaporation the part of active gluons do not convert into hadrons. They stay in QGS and 
become sources of soft photons (SP). Further we will analyze the experimental effect of 
SP excess (it is impossible to describe them by means QED). 
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Vs=91.4 GeV 






Figure 1: MD at 91.4 GeV. Figure 2: MD at 189 GeV. 



Figure 3: N g = aN q . 



10 20 30 40 50 



20 40 60 80 100 120 140 160 18t 



Figure 4: n h g = an q . 
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Figure 5: H q at 189 GeV. 
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Figure 6: cr(n) in pp. 



In the thermodynamic model without branches the active gluons which appear in the 
moment of the impact may leave QGS and fragment to hadrons. We consider that active 
gluons evaporated from QGS have Poisson MD with a mean multiplicity m. Using the 
idea of the convolution of two stages (fTTj) as well as the binomial distribution for hadrons 
from gluons we obtain MD of hadrons in pp-collisions in framework of the two stage 
thermodynamic model (TSTM): 



Pjs) 



m=0 



m! 



/ _/ l \rr 1 JV-( fl -2) 

1 -n) <"> 2 >- 



(19) 



P 2 ei (s) = e m . Our comparison (|19|) with the experimental data ^T] (see Fig. 8) gives 
values of parameters N = 4.24 ± 0.13, m = 2.48 ± 0.20, n h = 1.63 ± 0.12, and the 
normalized factor Q = 2 with % 2 /ndf ~ 1/2. We are restricted in sum ()19|) ME = 6 (the 
maximal possible number of evaporated gluons from QGS). The found gluon parameters 
iV and n h agree with the values of these parameters obtained at the e + e~ annihilation 
jH]. From TSTM the maximal possible number of charged particles is 26. This quantity 
is the product of maximal multiplicities of active gluons and of the maximal number of 
hadrons forming from one gluon ME ■ N. In TSMB there are no restrictions of this sort. 
Project " Thermalization" is partially supported by RFBR grant 03-02-16869. 
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